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During the last stage of collapse of a compact object into the horizon of events, the potential 
energy of its surface layer decreases to a negative value below all limits. The energy-conservation law 
requires an appearance of a positive-valued energy to balance the decrease. We derive the internal- 
state properties of the ideal gas situated in an extremely strong, ultra-relativistic gravitational field 
and suggest to apply our result to a compact object with the radius which is slightly larger than 
or equal to the Schwarzschild's gravitational radius. On the surface of the object, we find that the 
extreme attractivity of the gravity is accompanied with an extremely high internal, heat energy. 
This internal energy implies a correspondingly high pressure, the gradient of which has such a 
behavior that it can compete with the gravity. In a more detail, we find the equation of state in the 
case when the magnitude of the potential-type energy of constituting gas particles is much larger 
than their rest energy. This equation appears to be identical with the general-relativity condition 
of the equilibrium between the gravity and pressure gradient. The consequences of the identity are 
discussed. 
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Motto: " Chandrasckhar shows that a star of 
mass greater than a certain limit remains a 
perfect gas and can never cool down. The 
star has to go on radiating and radiating and 
contracting and contracting until, I suppose, 
it gets down to a few kilometers' radius when 
gravity becomes strong enough to hold the ra- 
diation and the star can at least find peace. I 
felt driven to the conclusion that this was al- 
most a reductio ad absurdum of the relativis- 
tic degeneracy formula. Various accidents 
may intervene to save the star, but I want 
more protection than that. I think that there 
should be a law of Nature to prevent the star 
from behaving in this absurd way." 

Sir Arthur Eddington [l[ 



I. INTRODUCTION 

After a star with a mass lower than the Oppenheimer- 
Volkoff limit has spent its nuclear fuel, it becomes either 
the white dwarf or neutron star. The rest of its ther- 
mal energy is irradiated into the neighbouring space and, 
thus, its temperature permanently decreases, down to a 
value approaching the absolute zero. Its internal struc- 
ture can be described with the help of theory of cold 
degenerated electron or neutron gas. 

In this work, we polemize with the current theory of 
the supplementary case: the end-state of the star with 
the final mass exceeding the Oppenheimer-Volkoff limit. 
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Namely, there are several serious arguments that such a 
star cannot end as a cold object. Even if it was cooled 
enough during the first stage of the collapse, its inter- 
nal energy and temperature can be expected to increase 
over all limits, when its radius approaches the horizon 
of events, which was found at the Schwarzschild grav- 
itational radius, R g . More specifically, the increase of 
internal energy can be expected when the proper radius 
of the object is reduced below, say, 1.01 R g or 1.001 R g . 

Our considerations are based on the classical quan- 
tum statistics and classical general relativity (with no- 
hair theorems valid), which both are well-known a long 
time. Only new element concerning the assumptions in 
this work is a term corresponding with a potential-type 
energy in the formulas of quantum-statistics equations 
of state for the extremely hot gas. In fact, our idea 
is not new, in principle. The potential-energy term is 
well known in the classical Maxwell-Boltzmann impulse- 
distribution law of ideal gases having a low internal en- 
ergy. In the latter, it is identical to the classical potential 
energy and situated in the Boltzmann part of the distri- 
bution law. 

In our description of gas properties, we use the 
formula for the relativistic energy which is suitable 
in a situation when the rest energy of gas particles 
can be neglected. In the classical works by Landau 
0, Chandrasekhar 0], or Oppenheimer & Volkoff [H, 
the influence of conservative gravitational field on the 
state properties of the gas was not taken into account. 
A step toward an inclusion of this influence occurred 
in more recent studies of stellar collapse and neutron 
stars. Prakash et al. @ considered the potential-energy 
term in course to provide the equation of state of a 
hot gas, which constitutes a neutron star a short time 
after its formation. The authors, however, considered 
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that approximative form of the relativistic formula for 
energy, which is applicable at the values much lower 
than the rest energy. Their work was focused especially 
on the nuclear reactions and alternative forms of matter 
constituting the star. They did not deal with any 
collapse going on down to the gravitational radius. Very 
recently, Moustakidis & Panos @, who continued in the 
work of Prakash et al., provided an equation of state 
for a hot nuclear matter with the /3-decay equilibrium. 
This kind of matter was again assumed to constitute 
a neutron star in less extreme conditions, allowing 
neutrons to be unstable, than we are going to consider. 

The following arguments can be presented to support 
the concept of the occurrence of extremely hot gas in an 
object with radius approaching R g . 

(1) Let us consider a radiation coming from an ultra- 
relativistic compact object to an observer at a relatively 
large distance. In his own reference frame, the observer 
detects the radiation which is red-shifted. In other words, 
its energy is lower in the reference frame of the observer 
than was at the moment of its emission, in the reference 
frame of the source. In a reverse situation, when a radia- 
tion comes onto the surface of the compact objects from 
a distant source, its wave-length is blue-shifted. This 
time, its energy is increasing with time. It is higher in 
the object-surface reference frame than in the frame of 
the distant source. 

The numerical calculation reveals that the frequency of 
the radiation coming to the horizon at R g from an outer 
space increases over all limits implying the corresponding 
energy to increase over all limits. Such the increase must, 
however, occur not only in the case of radiation, but in 
the case of non-zero rest-mass particles, when these are 
falling down, as well. The requirement of the local con- 
servation of energy in a volume of the compact-object- 
surface layer means that the increase must be, on the 
other hand-side, balanced by a decrease of that type en- 
ergy of the surface layer, which has a character of poten- 
tial energy. We refer to this energy as the "potential- type 
energy" . 

Therefore, when the collapsing-object radius ap- 
proaches the gravitational radius, its internal thermal 
energy must increase over all limits, because we do not 
know any efficient mechanism of the extreme-rate cool- 
ing. A fraction of internal energy is irradiated in form 
of photons and carried out by neutrinos only from the 
object's surface. Even in this case, the efficiency of both 
cooling processes decreases down to zero when the gravi- 
tational radius is approached (see our more specific con- 
siderations in Sect. IVlj ). Thus, an extreme accumulation 
of the internal energy seems to be inevitable in the last 
stage of the collapse to the horizon and can be the cause 
of a strong pressure with a steep gradient resisting to the 
gravity. After the object's radius reaches the horizon, the 
accumulated energy is trapped inside the object forever. 
(Some exceptions resulting from our work are discussed 
in Sect. IVIIl In the context of our work, we do not con- 



sider fine effects as, e.g., Hawking's evaporation of black 
holes.) 

(2) Tolman Q and Tolman & Ehrenfest [| found that 
inside the sphere filled in with an ideal material fluid, 
the product of both proper temperature of the fluid, T, 
and square root of component c/44 of metrique tensor is 
a finite constant, i.e. 

Ty/ga — Texp(^/2) = T C = const. (1) 

(exp(^) is an alternative expression of 344; see Eq.([2]) and 
the text below this equation) . It means that if a volume 
of the fluid is situated near the surface of the sphere, 
the radius of which approaches R g , then (744 — > and 
the temperature of the fluid must approach infinity the 
product of both could be the finite constant. 

(3) At the same time, Tolman & Ehrenfest Q demon- 
strated that the quantity P ra d 944 is equal to a finite 
constant inside the sphere filled in with the black-body 
radiation. P ra d is the radiation pressure. When a 
volume of gaseous material fluid sinks to a more and 
more curved space-time at the gravitational radius, 
radiation pressure P ra d must be magnified over all limits 
to keep the constancy of P ra d #44 • The requirement of 
this finite constancy also supports the hot concept of 
final stage of very massive stars. 

Considering all these facts, we search for the "hot" so- 
lution of the static configuration of a spherically symmet- 
ric object with no generated energy. Specifically, we give 
a brief summary of the theory and results that have been 
achieved up to date in Sect.|TTJ The introduction of a new 
modification of concerning fundamental theory with some 
obvious steps toward the solution of established funda- 
mental integrals is presented in Sect. IIIII In Scct. lIVl we 
derive the solution, which is applicable on an extremely 
hot gas. In more detail, we derive the state quantities of 
the ideal gas situated in a region of extremely curved 
space-time of the stable spherically symmetric object, 
whereby the curvature is so large that the magnitude 
of the potential-type energy of a given particle is much 
greater than its rest energy. In such a gas, a radiation 
must necessarily occur. So, an influence of the radiation 
pressure on the object's stability is discussed in Sect. fVl 
In Sect. I VII we attempt roughly sketch a behavior of the 
collapse of object before it disappears below horizon. Fi- 
nally, we summarize our most important results and dis- 
cuss some implications toward a potential observational 
evidence of this kind of objects, in Sect. IVIIl 

II. THE OLD THEORY 

The theory of the object contracting to approach a 
zero radius was completed by Oppenheimer & Volkoff 
They went on in the work started by Landau Q 
who showed that for a model of the star consisting of 
a cold degenerate Fermi-Dirac gas, after all the ele- 
ments inside available for thermonuclear reactions are 
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used up, there exist no stable equilibrium configurations 
for masses greater than a certain critical mass. In con- 
trast to Landau, who considered the Newtonian gravita- 
tional theory, Oppenheimer & Volkoff accounted for the 
general relativistic effects, which were expected for the 
very high masses and densities. 

Oppenheimer & Volkoff [H considered the static line 
element known for the spherical symmetry in the form 



ds 2 = 



-e x dr 2 - r 2 d9 2 



2 • 2 , 

r sin ( 



+ e"dt 2 , (2) 



where r, 9, and <fi are common spherical coordinates, 
t is time, and A, v are functions of r. They further 
considered the Einstein equations of field written for 
the case of spherical symmetry in which they put 
the expressions of the non-zero components of energy 
tensor T x : = T 2 = T 3 3 = -P gas and T 4 4 = E gas ex- 
pressed in this form by Tolman jlO| (see Eqs. (fTT)) . (fT5| . 
and (fl"9j). in Sect. 3). P gas and E gas are the gas pres- 
sure and its energy density (internal energy of gas per 
unit volume), respectively. The contra- variant compo- 
nents of the energy tensor for an ideal material fluid, 
T n = T 22 = T 33 = p and T 44 = were specified 

earlier by Eddington Q. 

To solve the field equations, Oppenheimer & Volkoff 
established the auxiliary function u = u(r) equal to 



-r(l 



(3) 



with the help of which the equations of field could be 
written in form (in SI units used throughout our paper) 



du 
dr 



AttG 



-E, 



dP n , 



E„ 



dr 



r(r — 2u) 



gas' 



PgasT 3 + U 



(4) 



(5) 



where G is the gravitational constant and c is the ve- 
locity of light. For P gas and E gas , they supplied the 
parametric form of the equation of state, for the cold de- 
generated Fermi-Dirac gas, found by Chandrasekhar (0], 
appendix) . After the analysis of the possible solutions of 
final equations, Oppenheimer & Volkoff concluded that 
there was not any stable solution for the object having 
the mass exceeding the finite limit they specified. 

If their conclusion was indisputable, our re-opening the 
problem on the stable configuration would be meaning- 
less. However, we find three of their assumptions dis- 
putable. At first, they assumed that the collapsing ob- 
ject is cold enough for the theory of the cold degenerated 
Fermi-Dirac-statistics gas to be applicable. As we argued 
in Sect. HI the decreasing potential- type energy should be 
balanced by an increase of the thermal energy and this 
process leads, instead, to a hot gas. The extreme increase 
occurs especially at the last stage of collapse in proper 
time, but above the event horizon, when the physical 
radius of collapsing object is close to R g and metrique- 
tensor components g\\ and (744 start significantly deviate 



from unity. So, even if the object was cold enough at a 
proper radius several times larger than R g , it must be 
highly heated during this last stage. 

Interestingly, Oppenheimer & Volkoff noted, in their 
paper that for those singular solutions in which 344 
vanishes, it is conceivable that the temperature may be 
high with the respect to the Tolman & Ehrenfest Q rela- 
tion Tyjgn = const. And, the equation of state reduces 
to E gas — iP g as in the limit of T — > 00 (which we anal- 
yse in Sect. II VI) . However, they supposed the appropri- 
ate conditions for the very high temperature only in the 
center of the object. They did not noticed that the ap- 
propriate conditions may also exist in the region around 
the border of the sphere of radius equal to R g . 

At second, they required the pressure approached zero 
at the surface of the object, i.e. for r = Rb (Rb is the 
proper radius of the object, i.e. the distance of the ob- 
ject's surface from its center measured in the object's 
surface reference frame). In every-day experience, we 
know, however, a lot of examples of pressure or density 
discontinuities: the border between a rocky earth surface 
and air, the sea level and air, etc. Another example can 
be found in the outer space, where neutron stars are ex- 
pected to have a very large density at the surface, but, of 
course, zero density above the surface. We cannot apriori 
exclude the possibility that the pressure (or energy den- 
sity or particle concentration) is finite in the surface layer 
between radii Rb — dr and Rb and zero in the adjacent 
layer above, between Rb and Rb + dr. 

At third, they started the numerical integration of the 
final equations with the initial value u Q — u(r — 0) being 
in the interval > u Q > —00. They excluded the pos- 
sibility u > arguing that the component of metrique 
tensor gn — — e A must never be positive, therefore u Q = 



for all finite values of e 
ues of e _A at the origin. 



and u a < for infinite val- 
The argument of not positive 
is, however, questionable. Let us consider the ob- 
ject which already collapsed below the gravitational ra- 
dius. In the interval of distances between Rb and R g , 
the well-known outer Schwarzschild solution is applica- 
ble. According to this solution — e A = —1/(1 — Rg/r) > 
for Rb < r < R g . Therefore, if we rejected the possibil- 
ity — e A > 0, there would not exist not only any solution 
for the stable configuration of an extinct star, but we 
would not have any acceptable description of the field 
for r < R g , either. Consequently, the reason for the in- 
finite collapse due to the gravitational attraction would 
disappear. But we believe that it is reasonable to retain 
the possibility — e A > 0. And, if we accept this possibility 



for r < R a 



then why we should not doing so for r > R g ? 



III. THE BASIC CONSIDERATIONS 

Let us consider a material object, which is spherically 
symmetric and ultra-compact, i.e. its proper physical ra- 
dius, Rb, is only slightly larger than its Schwarzschild 
gravitational radius, R g . To better explain the conse- 
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quences resulting from the new formulas we find below, 
we consider the idealized, non-rotating object that is in 
or near its supposed stable configuration. Hereinafter, 
we refer to this object as the " ultra- relativistic material 
sphere" (URMS). 

In this section, we summarize the well-known funda- 
mental equations of the problem of URMS stable con- 
figuration and describe several first steps to obtain the 
solution for the stable URMS. Within the summary, we 
introduce our own modification of the fundamental equa- 
tions necessary for taking the strong, ultra-relativistic 
gravitational field into account. The steps toward the 
solution can be made generally for the gas consisting of 
a single kind of particles having spin 1/2 (e.g. neutrons). 
No other assumptions are needed. 

In the following, we consider the reference frame re- 
lated to a volume of URMS' surface layer. If the URMS 
is stable, then the number density, n gas , pressure, P gas , 
and internal energy per unit volume, E gas , of the gas 
in the surface layer can generally be calculated, as the 
functions of the impulse p, with the help of well-known 
integrals (e.g. (Tlj |) 

n gas = / n gas {p)^p 2 dp, (6) 



Pgas = - / n gas (p)pVAiTp 2 dp, (7) 



Egas = I n ga s{p)£(p)4:TTp 2 dp, (8) 
JP 

where the kinetic energy £ and velocity V of a particle 
can be calculated by the well-know relativistic formula 

£{p) = Vc 2 P 2 + W 2 - W , (9) 



Wj/(kT) of the argument in Eq. (fTTj) must converge to ex- 
pression p 2 /(2mkT) + W p / (kT) figuring in the Maxwell- 
Boltzmann distribution law. So, the energy does not 
contain the rest-energy term, W . It is only the sum 
of kinetic and potential-type energy, i.e. 



W k+P = ^c 2 p 2 + W 2 -W ~x, (12) 

where — x is the potential-type energy of the particle. By 
the convention used in our work, the fact that the gravity 
is attractive force is expressed with sign minus in front of 
X, in Eq. (|12p . Within this convention, it is always valid 
that x > 0. Completing the description of the energy 
state with term —\ is an essential point of our work. In 
Sects. HVHVIll we derive some important consequences 
of this additional term. 

The relativistic formula (fT2")> . but with the rest-energy 
term included, is well-known in quantum physics. Its va- 
lidity was proved, for example, in the construction of the 
Klein-Gordon equation, which is the basis to derive the 
well-known Dirac equations providing the energy terms of 
hydrogen atom, which are in a very good agreement with 
the experimental values. In this usage, — x is identified to 
the classical Coulomb energy for a static system. In the 
strong, ultra-relativistic gravitational field, the function 
X — x{ r ) cannot be, of course, identified with the well- 
known formula for the potential energy in the Newtonian 
gravitational field. The explicit form of this function for 
a strong field is unknown. We discuss its dependence on 
the quantity v (in metrique tensor; see Eq.([2])) in the case 
of the extremely strong field in the next section. For now, 
we profit from working with the implicit x- Fortunately, 
a more specific expression of x appears not necessary. 

In the context of function Xi the- "Fermi impulse", pp, 
corresponding to the Fermi energy, is 

PF = -V(x + W + f i) 2 -W 2 . (13) 



V 



d£ _ c 2 p 

dp ^c 2 p 2 + W 2 



(10) 



Symbol W stands for the particle's rest energy. 
The number density n gas {p) for the ideal gas is 



fi 



(p)-n]/(kT)} + V 



(11) 



where Kn is a normalization constant, k is the Boltz- 
mann constant, T is the temperature of the gas, fj is the 
degeneracy of state j (i.e., it gives the number of states 
having the same energy Wj), and /j, is the chemical po- 
tential. 

The energy state of a particle Wj in the argument 
of exponential in Eq. ljlip must be of the type varying 
with impulse p and with the force field in which the 
gas is situated. In the limit of a weak field, the part 



This relation can be derived requiring the argument of 
exponential in Eq. (|ll[> to be zero and identifying the en- 
ergy Wk+p given by Eq. (fl"2|) to the Fermi energy of the 
particle. We note that pf — x/ c f° r X ^ W and, of 
course, pf — * oo for x ~> °°- 

The energy states, Wj, in Eq. lfTTj) differ each other by 
their rest energies, W j, for various kinds of the particles. 
Since we consider, in this section, only the gas consisting 
of one-half spin particles of one kind (neutrons), the sum 
in Eq. ljlip consists of only a singe term (j = 1) and fi = 
2. 

The size of impulse p can vary from zero to infinity. 
So, the number density, pressure, and internal energy 
per unit volume are the integrals 



' go s 



= K 



P 



exp 



y / c 2 p 2 +W^-W - X -l-t 

FT 



-dp, (14) 



1 
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p — 

r gas — 



Kc 2 



o ^/cW+W 2 



exp 



\/cV+W -w -x-v 



kT 



dp, (15) 



E q as = K 



p 2 (yjc 2 p 2 + W 2 - W ) dp 

( Jc 2 p 2 + W 2 -W - x -ti 
ex P FT 



'SPgas - W n gas + KW 2 



o VcV + w 2 



exp 



vW+W - W - x - /i 



dp, (16) 



respectively. We denoted K — 8irKjy. 

For the URMS, the Einstein's equations of field can be 
written in their form for the spherical symmetry Q , [Ic| 
(see also [HI, §95): 



— P = exp(-A)(-- + - 



— r P = exp(-A) 



(17) 



1 dX dv 
4 dr dr 
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E = exp(-A) 



1 dX 

r dr 



(18) 



(19) 



in the SI units. P and E are the total pressure and 
total energy density in a given volume of space in the 
proper reference frame. The cosmological constant, A, is 
put zero, since we assume that the URMS can occupy 
a volume many orders of magnitude smaller than the 
cosmological distance-scale. 

It is well-known [l(| , [13] that we can derive the equa- 
tion 



is the mass within radius r. There is an essential differ- 
ence between the relativistic Eq. (|20|) and classical-physics 
Eq. (|2T|) . The latter does not contain any term corre- 
sponding to term —{P/2)dv/dr in Eq. (l2"Ul) . Recently, 
Rappaport et al. [1 31 ] and Schwab et al. [14| represented 
this term as the self-gravity of pressure. 

Our task to find a stable configuration of the URMS 
will be solved, when we find a reasonable combination 
of functions P, E, and v, for which Eq. ([20f is satisfied. 
For the gas, we calculate the corresponding, gas-related 
part of the left-hand side of Eq.([2"0]) derivating Ea.lfTS]), 
in which the temperature, T, and potential-type energy, 
-X, are expected to be the functions of radial distance, 
r, from the center of the URMS. (Since we assume the 
mono-particle gas, the chemical potential, /i, is obviously 
not any function of r.) We obtain 



UJ gas 

dr 



Kc 2 dT 
ikT 2 dr 

Wo 



Ipi 



Kc 2 
3kT 
V[idT 



dr 



d_X_ 
dr 

IP2, 



(22) 



where 



'pi 



I Jc 2 p 2 + W 2 -W - X -H 

oo P exp ( w 



Jc 2 p 2 +W 2 -W -x-» 

exp ( * w 



-jdp (23) 



and 



<P2 



4 / y/c 2 p 2 +W 2 -W - X -V 

p exp ( * w 



y/c 2 p 2 + W 2 



^c 2 p 2 + W 2 -W -x-A 
kT j 




/ ^Jc 2 p 2 + W 2 -W -x-H , . , 

ex P I FT I + 1 



dp. (24) 



dp, (25) 



dP 

dr 



P dv 



dr ' 



(20) 



from Eqs. ([T7| — (fT9| . This is the relativistic analogue of 
the Newtonian relation 



dP 
dr 



dr 



GM r 



(21) 



giving the dependence of total-pressure gradient, dP/dr, 
on the gravitational potential, \I/ ([12], §95). This equa- 
tion is also known as the equation of the hydrostatic 
equilibrium in the theory of internal structure of stars 
(e.g. [ill, chapt. 1). Quantity p is the density of stellar 
plasma at the distance r from the stellar center and M r 



S F 



P 



o VcV + W 2 
-i -l 



and 



/ y/c 2 p 2 + W 2 ~Wp^X-A a I 
CXP kf +1 



V 



dp, (26) 



o ^c 2 p 2 + W 2 
-, -l 



/ y c y + w 2 - Wo - x - M , , 
exp kf +1 



dp, (27) 
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integrals ([23]) and (|24[) can be expressed (after one step 
of the per partes integration) as 



Jpi = AkTSp 



3kT 

Ip2 — — 5~ S n - 



(28) 



(29) 



Since the integrals S n , Sp, and Se can also be given as 
S n = nga S /K, Sp = 3Pgas/(Kc 2 ), and Se = (E gas - 
3P ffas + W n gas )/(KW*), the gradient of pressure 
can be expressed as 



potential-type energy, \- Consequently, the relation ([12]) 
reduces to 



W = cp — x- 



(33) 



As well, we neglect the chemical potential, /i, because any 
chemistry cannot be important in the considered extreme 
conditions. 

If both W and \x are neglected, integrals 1)14 — (|16 p 
are simplified to 



K 



P 



exp(^) + l 



dp, 



(34) 



dPna 



gas 



dr gas dr T dr 

In the last equation, we denoted 



P* Egas ~t~ Pgas X^gas ^^gas 



(30) 



(31) 



IV. THE SOLUTION FOR AN EXTREMELY 
HOT GAS 

Let us consider the reference frame, to which the line 
element ([2]) is referred. It means, it is the proper ref- 
erence frame of an infinitesimally small gas volume in 
the URMS surface layer at the distance r from the cen- 
ter of URMS. We further consider the static case, i.e. 
the components of the metrique tensor gn = exp(A) and 
#44 = exp(^) are the functions of only the space coor- 
dinate r. According the outer Schwarzschild solution of 
the Einstein's equations of field, it is well-known that 
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1 

344 



1 



1 - Rg/r 



For 



R g , there appears the singularity \l/g. 
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00, which is essential in our deductions. The behavior 
of the state quantities of gas are derived in the coor- 
dinate system, in which this singularity appears. Its 
existence implies that the conditions for the extremely 
high pressure and temperature must exist not only at 
the center of the URMS (at r — > 0, therefore 1/r — ► 00), 
but around the sphere of radius R g as well. This indi- 
cates that it is reasonable to investigate, in the search 
for the URMS' stability, the properties of the extremely 
hot, ultra-relativistic gas, the existence of which appears 
to be very probable in the surface layer of the URMS 
with the proper physical radius not much different from 
its gravitational radius. 

Close to R g , we can obviously expect that x ~~ * 00 ■ So, 
let us search for the solution of integrals (fbl]) — (TIB)) as- 
suming that the gas is situated in the gravitational field 
characterized with the strong inequality x ^ W Q - In 
other words, we neglect the rest energy of the gas par- 
ticles with respect to the magnitude of their relativistic 



Pa. 



Kc 
E„ 



exp 



dp, 



3P 

lJ± gas ) 



(35) 



(36) 



respectively. One can notice t hat we also neglect W in 
the argument of square root <J c 2 p 2 + W£ in the denomi- 
nator of Eq. (|15p , where the rest energy W Q is comparable 
to cp at the down limit of this integration. Fortunately, 
this neglection is possible, because the contribution of 
the integration in this region of p to the total integral is 
well-negligible. 

Eqs. lfM)) and ()55)) do not contain any particle charac- 
teristics, except of constant of proportionality K. Since 
all following relations arc valid for every constant K, our 
considerations concerning the hot gas can be generalized 
for an arbitrary Fermi-Dirac-statistics gas consisting of 
several kinds of particles characterized with different de- 
generacy states fj and having the partial concentrations 
Ug t j. In the generalization, the constant K = 47r/i, for 
(32) the 1/2-spin mono-particle gas, must only be replaced 



with constant K — (47r/n 9as ) 
half-spin multi-particle gas. 
The integrals of type 



fj for the various- 



7 T = 



exp(a; — b) + 1 



-dx 



(37) 



are well-known in quantum statistics. The approximative 
solution can be found in form (e.g. (l5j ) 



3 =J ^n( 7 - 2j + 1) 
n( 7 ) 



11(7 + 1) 



n(«7), 



(38) 



where 



f (-1)' 

1=0 y ' 
\TZ(b,-y)\ <exp(-6), 
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and II is the factorial function. Below in this section 
we however find the possibility that exponent b in the 
condition delimiting \TZ\ may not be large, therefore \1Z\ 
cannot be neglected. For the values of 7 figuring in the 
considered integrals (f34| and (|35]) . i.e. 7 = 2 and 3, we 
utilize the full results of integrations yielding 



l gas 



K X 3 
3c 3 



■E 



12a 2 (kT) 2 6(kTf 



.r 



cxp —j 



— ) 

kTJ 



(39) 



9as 12c 3 



24a 2 (fcT) 2 48a 4 (fcT) 4 



24(fcT) 4 



E 



exp 



, (40) 



where a 2 = tt 2 /12 = 0.822467033 and a 4 = 0.947032829. 

Using Eq. (136|) and neglecting fi, the auxiliary quantity 
B (see Eq.^P) is 



^Pgas X^gas 

Using further Eq.fl]), we can find 



1 dT 



1 dv 
2~dr' 



With the help of Eqs.gTj) and 
relation (l30l) in form 



(41) 
(42) 

we can write the 



dPna 



dv 



dr 



d\ 1 dv 
dr 2^ dr 



(43) 



For the strong field characterized with condition \ ^ 
W Ql function \ — x{ v ) can be found assuming that the 
URMS was formed as a product of collapse of an object 
with so large radius, that the magnitude of potential- 
type energy of gas particles was much lower that their 
rest energy, i.e. x ^ Wo for a single particle. And, it 
was not extremely hot, therefore the kinetic energy of the 
particles was also much lower than their rest energy, i.e. 
yj c 2 p 2 — W 2 — W <C W for a single particle. It means, 
whole energy of a particle consisted only of the rest en- 
ergy and we can, in the given context, well approximate 
the sum of kinetic and potential-type energies by zero 
(Wk+p = 0) before the collapse. The kinetic energy per 
one particle can also be given as E gas /n gas . Since the 
total energy is expected to be conserved, energy Wk+ P 
must be zero also after the collapse in the situation, when 
X W . This requires the validity of E gas /n gas — \ = 0. 

The energy density, E gas , is proportional to the tem- 
perature. Specifically, we express it as 



E„ 



sn 



•AT, 



(44) 



where s is a constant parameter. Using relation (|44|) and 
fact that E gas /n gas — x = 0, function x is 



X = skT. 



(45) 



If we now use the Tolman & Ehrenfest's relation HJ, we 
obtain 



X = skT c exp(-v/2). 



(46) 



The constant skT c can be identified to W because of 
the following requirement. Let us consider the quantum- 
physics concept of the wave associated with the particle. 
If the angular frequency of this associated wave in the 
free space, far from the URMS, is uj , the wave-related 
energy of the particle is identical to its rest energy and 
given by the well-known formula W a = huj , where h is 
the Planck constant divided by 2ir. In the place of a 
considered gas volume, where the potential-type energy 
is x, the angular frequency is lu and the corresponding 
wave-related energy of the particle is W w = hw. Energy 
W w is positive and increases with the increasing magni- 
tude of the gravity. It is obviously a consequence of the 
fact that just this increase of energy balances the corre- 
sponding decrease of the particle's potential-type energy. 
Therefore, we can put x = W w — W . Using further 
the well-known transformation between w and u> , i.e. 
lu = io j ^/g44 , we finally find 

W w - W = x = exp(-z//2) - huj , (47) 

or, after neglecting W Q , 



X = Wo exp(-i//2), 



(48) 



where we put ftw — W Q . Comparing the right-hand side 
of this equation with the right-hand side of Eq. (j46|) . it is 
clear that constant skT c must equal W . 

In the last step in the search for the URMS stabil- 
ity, let us look for such a configuration of quantities that 
the quantum-statistics solution for the gradient of pres- 
sure satisfies the equation of the equilibrium (|2"D|) . When 
Eq. (|36|) is also the valid relation between the total inter- 
nal energy density and total pressure, i.e. it is valid that 
E = 3P (this validity can be deduced from Eqsl3"6"]) and 
([Ml)), Eq.® changes to 



dP = _ 2p dv 
dr dr 



(49) 



We can see that this functional form for the gradient of 
pressure, derived from the Einstein's field equations, is 
identical to that derived within the quantum Fermi-Dirac 
statistics (see Eq. (l4*3^ ). if 



dx 
dr 



1 dv 
2 X dr~ 



0. 



(50) 



The validity of this condition can easily be shown differ- 
entiating Eq. (|4"8"l) . therefore Eq. (|4"3")) becomes simplified 
to 



dP ga 
dr 



-IP 



dv 
' dr 



(51) 



8 



If the radiation was completely absent, this would prove 
that the URMS would be close to its stable configuration 
at whatever value of the component exp(V) of metrique 
tensor or whatever radius Rb for which the condition 
X ^> W (or W /x <C 1) is satisfied. (The exact sta- 
bility would appear for W /x ~~ * 0. A further discussion 
of the stability is given in Sect. I VII ) 

The continuity of the gravity behavior implies that the 
function exp(— v) as well as the derivative du/dr at the 
surface of the URMS should approach their values cal- 
culated by the outer Schwarzschild solution, according 
which 



exp(^) = 1 



Ra 



and, subsequently, 



dr 



(52) 



(53) 



We can see that du/dr — > oo for r 
dPgas/dr — > oo (see Eq. (|5Tj) for r — ► R g 
Integrating Eq. (f5Tj) . we can find 

Pgas = P exp(-2v), 



R g , therefore also 



(54) 



where P is an integration constant. This solution trans- 
parently implies that P ga s —> oo when exp(— v) — > oo. 
Constant P a can be found realizing the fact that the so- 
lutions (|4"0|) and ([M]) must be identical. Analysing these 
functional forms, we can find 



Po = 



24 
s 



KW |\ , 24a 2 , 48a 4 

4 2^ ~l eX P(--? S ) 

3=1 



3 



where 



kT r . 



const. 



(55) 



(56) 



with respect to the relations fl} and (|48p. 

For the cold, degenerated relativistic Fermi-Dirac gas, 
the relation between the pressure, P gas , and concentra- 
tion of particles, n gas , was found in the form P gas oc n^as 
(see, e.g., [IH). It is interesting that this proportionality 
is also valid in the case of the extremely hot gas. From 
Eas.pg]) and we can obtain 



we can re- write Eg. (139)) as 



'"gas W>oX ? 



(59) 



from which x = ("■gas/^o) 1 ^ 3 - Supplying this into 
Eq. ([5*7|) . we obtain 



p — 

r gas — 



Pn 



,4/3 = 



4/3 9 as 



const, x n 



4/3 
gas- 



(60) 



V. THE INFLUENCE OF RADIATION 

Since the strong gravitational field of the URMS im- 
plies a high internal energy of gas and, therefore, a high 
temperature, we can expect an occurrence of an intensive 
radiation inside such a gas. The pressure of this radiation 
has obviously to contribute to the total pressure. Let us 
estimate an influence of this contribution on the stability 
of the URMS, in this section. 

The statistical behavior of photons of electromagnetic 
radiation can be described by the Bose-Einstein statis- 
tics. Because of a larger context, we give a review of the 
well-known derivation of basic relations for an arbitrary 
kind of bosons. 

In an analogy to Eqs. ([34|) — (|36|) of the Fermi-Dirac 
statistics, the concentration, pressure, and internal en- 
ergy per unit of volume can be written as 



riBE = K 



exp(^)-l 



BE 



Kc 



exp 



, kT > 



dp, 



-dp, 



E BE = Kc 



o exp 



' cp-x ' 
v kT I 



- 1 



dp = 3P B e, 



(61) 



(62) 



(63) 



for the photons, or other Bose-Einstein-statistics parti- 
cles, when their rest energy, W , is zero or can be ne- 
glected. Again, we inserted the potential-type-energy 
term, — x, into the arguments of the exponentials to com- 
pletely characterize the energy state. 

Calculating integrals (|6"Tj) and (|62"j) . we can find 



tibe 



2K ^ i 



3=1 J 



exp(i^), (64) 



Denoting further 
K 



Pn, 



12a 2 



( _l)3+i 



s 3 

3=1 



(57) 



•exp(-js)] 



(58) 



Pbe 



C 3=1 3 



kT) 



(65) 



The solutions ([64)1 and (|65|) for nsE and Pbe contain 
the sums, which are divergent for xi ' {kT) > 0. In an ex- 
pectation of a stable gas volume, this result implies that 
either the same initial assumptions cannot be valid for 
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both kinds of particle statistics, Fermi-Dirac and Bose- 
Einstein, or x/{kT) = 0. It is actually known that the 
bosons are not permanent constituents of matter. In ex- 
periments, they are found to be unstable, except of pho- 
tons, therefore no law of the conservation of their number 
is satisfied, in a common reality. 

For photons with zero rest mass and, thus, zero rest 
energy, Eqs. (|61|) — (|63|) are exact, not only approximative 
equations. Consequently, the result of integrations ex- 
pressed by relations and (|rJ5|) is exact. Eq.([i8"|) im- 
plies that x — for the photons, since W Q = for them. 
As mentioned above, the sums in relations and (|6"5)) 
are convergent for x — 0, therefore we have the regular 
solutions for Ube, Pbe, as well as Ebe (Ebe = 3Pbe)- 

For the radiation, the normalization constant Kpj is 
equal to Kn = l/(2irh) 3 . The latter implies that 
K = l/(TT 2 h 3 ). If x — 0, the pressure given by 
Eq.§5\) becomes P BE = Prad = 2(3 4 (kT) 4 /{ir 2 h 3 c 3 ) = 
Tr 2 (kT) 4 /(45h 3 c 3 ) (where [3 4 = E^ii -4 = tt 4 /90). 
Or, with the help of the radiation density constant, 
a = 7r 2 fc 4 / '(15h 3 c 3 ), we can rewrite this relation to the 
form of another well-known formula for the radiation 
pressure 



Prad = 7;aT 4 . 



The gradient of this pressure is 

dPrnA B r „d dT 



where 



dr T dr 1 



Brad — 3 fl ^ 4 ' 



(66) 



(67) 



(68) 



Using Eqs.(fr]) and (f6"6")h the quantum statistics Eq. (|6"7| 
can be re-written to 



dPrad 

dr 



dv 



(69) 



i.e. to the functional form identical to Eq. ((49|) derived 
from the Einstein's equations of field. As well, it is iden- 
tical to the functional form in Eq. (|5Tj) valid for the gas 
consisting of fermions. We see that the gradient of the ra- 
diation pressure behaves in the same way as the gradient 
of the fcrmion-gas pressure. 



In terms of mathematics, the equation of the "asymp- 
totic" equilibrium between the total pressure gradient 
and gravity is the sum of Eqs.([5Tjl and (f69"|) . i.e. 



d(Pgas ~t~ Prad) 

dr 



2(Pgas 



dv 

Prad ) ~t ■ 

dr 



(70) 



We emphasize that this-type requirement of the stable 
configuration of URMS is set by Eq. (T49"|) yielding from 
the Einstein's equations of field. Our proof of the bal- 
ance (in the limit W /x — ► 0) is the fact that the same 
equation as Eq. (|49|) also occurs as the necessary descrip- 
tion of state quantities in the quantum statistics, which 
is represented by Eq. ([70"|) . (We note P gas + Prad = P, in 
this equation.) Both gravity and internal quantum state 
of the considered gas/radiation converge to the same be- 
havior close to the horizon, at the gravitational radius. 

When applying the above conclusion about the bal- 
ance in reality, we must, however, keep in the mind that 
the condition of the balance is x ^ Wo, or W /x -C 1. 
The exact balance is set, when W /x ^ 0. This can 
be expected for r — > R g . For a larger r, we can expect 
a minute break of the exact balance: in contrast to the 
before-URMS-stage fast collapse (according to the the- 
ory of the last stages of stellar evolution), the URMS 
either very slowly expands or collapses. If an expansion 
appeared, it could be fuelled only by the internal energy 
of gas. Since there is still a certain energy loss (see our 
considerations below), such the expansion would earlier 
or later turn into a the collapse. 

In this section, we discuss some aspects of the last stage 
of the assumed final, slow collapse of the URMS. We re- 
gard as the last that stage of the collapse, when quan- 
tity exp(— v) 1, but before the moment when all the 
URMS' matter disappears under the horizon of events. 
Unfortunately, we do not know the exact solution of gen- 
eral equations (|T3|) . (Tl5|) , and (fl7j|) giving the quantities 
of state as the functions of distance. Moreover, we would 
need a solution of Einstein's equations of field for non- 
static case, which is also unknown. A certain information 
can be deduced from a rough comparison of relative ratios 
of the internal-energy loss of two URMSs with different 
masses due to the escape of the energy to the outer space. 

Considering the relation E ra d = 3P ro d as well as 
Eq. (|6"6")) , in which T — T c exp— v/2 with respect to 
Eq.([T]), the energy of radiation in a unit surface volume, 
U 1; is 



VI. THE LAST STAGE OF THE COLLAPSE 
ABOVE THE EVENT HORIZON 

Since the individual gradients of both gas pressure and 
radiation pressure acquire the values approaching those 
which are needed to balance the gravity, as we demon- 
strated in Sects. IIVI and [V] the sum of these gradients 
also approaches that at which these gradients would bal- 
ance the gravitational attraction in the extremely strong 
potential close to the Schwarzschild gravitational radius. 



IU rQd = y 1 aT c 4 exp(-2^). 



(71) 



The energy irradiated through the unit surface area, Si, 
of this volume per a time unit can be given with help of 
the well-known radiation law as 



$ = SifiT 4 = SioT 4 exp(-2^) 



(72) 



where a is the Stefan-Boltzmann constant. 

When we assume that the URMS' gas consists of neu- 
trons with the rest mass m n , we can calculate T c using 
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Eq. (|56|) . i.e. T c = m n c 2 /(sk). In the thermodynami- 
cal equilibrium, the temperature of gas and radiation is 
the same. In our ultra-relativistic case, s = 3, therefore 
T c fa 3.63 x 10 12 K. 

Eq. (j?2"|) gives the radiation flux in the coordinate frame 
of the emission source, i.e. in the frame of a volume in 
the URMS surface layer. A distant observer detects the 
original radiation, having the angular frequency ui, with 
the red-shifted angular frequency lo , whereby 



lo = lo exp(z;/2). 
So, the flux measured by the distant observer is 



SioT c exp \ --u 



(73) 



(74) 



The fraction f rac [ of radiation energy detected by the 
distant observer per time t, which can be regarded as the 
loss of internal energy of the surface volume V± per time 
t. is 



frad = 77— texp(z//2), 

Via 



(75) 



in the simple case when t is short enough that f rac i -C 1. 
To provide a specific numerical estimate of f ra d, we would 
need to specify volume V\ . Unfortunately, it is impossible 
without an explicit knowledge of functions A = A(r) and 
v = v(r) for the URMS' interior. Nevertheless, we can 
expect the continuous behavior of quantity v. At the 
URMS surface, it must approach to that known from the 
outer Schwarzschild solution, i.e. exp(z/) = 1 — R g /Rb- 
Therefore, exp(f/2) — > for Rb — > R g . For a finite 
time t, the fraction of the radiation-energy loss must also 
frad — > 0, when R — ► R g . The URMS collapse can be 
expected to be stopped at Rb — R g . 

In a weak gravitational field with no significant red- 
shift, both internal energy per a unit volume and irra- 
diated energy from the volume per a unit time are pro- 
portional to T 4 , therefore the appropriate f ra d does not 
depend on temperature. There is no such the decrease 
in the efficiency of cooling as in the above mentioned, 
exp(^/2)-proportional, ultra-relativistic case. 

We can expect a certain dependence of f ra d and, thus, 
the efficiency of cooling on the mass of the URMS. To 
compare, at least very roughly, this dependence for the 
URMSs of two different masses, let us consider a naive 
concept based on the assumption that the whole internal 
energy of URMS is uniformly distributed in the volume 
A-kRI/Z. Considering further the corresponding surface 
47ri? 2 and identifying Rb = R g (the difference between 
Rb and R g is a small fraction of R g during the whole last 
stage of the URMS collapse), we can put S1/V1 ~ 3/R g 
implying f ra d ~ 3atexp(u/2)/(aR g ). Or, the time-scale, 
on which fraction f rac i of radiation energy is lost, can be 
estimated as 



O'Rg frad 

3ct 



exp(-i//2) « ^2lMexp(-i//2). (76) 



The whole context of our study intuitively implies 
a replacement of the old concept of the black hole by 
the concept of the URMS. Thus, not only stellar-sized 
objects, but also the massive compact objects in the 
galactic centers can be identified with the URMSs. The 
linear dependence of the time-scale t on mass M given 
by Eq. (j76|) indicates that the time-scale is about 10 5 
to 10 10 times longer for a super-massive URMS in a 
galactic center, with correspondingly larger mass, than 
that for a stellar-sized URMS. This large factor indicates 
that the collapse of the potential large URMSs in the 
universe cannot obviously be the matter of a moment or 
a short period. 

Another potential mechanism of the internal energy 
loss is the escape of neutrinos carrying energy. This 
mechanism is well-known to occur at the first stage of 
the collapse of a massive stellar object. It is verified 
observationally, so far, by the detection of neutrinos pre- 
ceding the optical outburst of supernova 1987A in the 
Large Magellanic Cloud [TtI ]. 

Although the extremely high pressure inside the 
URMS is keeping the particle constituents of the gas, 
neutrons, stable, the neutrinos can appear as a final re- 
sult of collisions of photons, which are abundant and very 
energetic in the URMS environment. Specifically, the 
collision of two photons can produce an electron-positron 
pair, which can again annihilate to produce the ncutrino- 
antincutrino pair. Or, neutrino (antineutrino) can occur 
in a collision of photon with a lepton (occurred in the 
first reaction), which behaves as a catalysator. Another 
reactions, including other elementary particles, can also 
occur in a collision of photons. The important fact in 
our context are the photons entering the reaction and 
neutrinos/antineutrinos as the result of the reaction. 

The mean free path of neutrino in the gas with con- 
centration of screening particles n gas can be given as 



1 



^u^gas 



(77) 



If the energy of neutrino is £„, the cross-section in the 
collision of a neutrino with a particle of the gas can be 
estimated as [llf 



t f 



(78) 



where the constant of proportionality, in SI units, was 
found to equal a ov ~ 10~ n m 2 J -2 . 

The neutrino cooling can only be sufficient if the energy 
carried out by a single neutrino is comparable or higher 
than the internal energy of gas per one gas particle (neu- 
tron), i.e. E V ^E '/n gas = skT. Supplying Eqs. (f59|) and 
(|78p into Eq. (|77p . whereby \ is expressed using Eq. 
we easily find for the neutron gas 



I// . 



1 



<J vn W% 



exp(5^/2) - 10- lb exp(5iV2) (79) 



(in meters). Identifying constant K to K = l/(ir 2 h 3 ) in 
relation (1581). the value of constant n Q is n = 2.253 x 



11 



10 75 m -3 J -3 . Since exp(5z^/2) <C 1, the mean free path 
of neutrino l v -C 10~ 16 m. This path is much shorter 
than the scale Rb — R g of the extreme-condition region, 
in which the reciprocal metrique tensor 1/(744 1. The 
neutrinos with a non-negligible energy can escape only 
from the upper-most, extremely thin surface layer to- 
gether with the photons. Since they can originate only 
from photons, an outflux of neutrinos cannot exceed the 
outflux of photons. The appearance of neutrinos (by con- 
version from photons) only means that a part energy, 
otherwise irradiated by photons, is carried out by neutri- 
nos. 

The efficiency of the neutrino cooling of the thin sur- 
face layer correlates with that of radiation. The neutri- 
nos also loose their energy during their escape from the 
URMS surface due to the red-shift as photons. And, ei- 
ther these particles cannot escape when the URMS radius 
reaches the gravitational radius. 

VII. SUMMARY AND CONCLUSION 

In Sects. ITVl and IVl we found the dependence of state 
quantities on the auxiliary function v, for an extremely 
energetic environment. This dependence on v appears to 
be the same for both particles satisfying the Fermi-Dirac 
statistics and photons satisfying the Bose-Einstein statis- 
tics. In the case of the high energies, the rest mass of 
fermions is negligible, therefore they seem to be as mass- 
less particles as photons, and the identical dependence 
on v is not very surprising. 

In the beginning of Sect. IVI1 we demonstrated that 
the behavior of the gradient of total pressure, given 
by Eg. ([70]) . is identical to the relativistic equation of 
state (|4l?)) . being a simplification of Eq.lJ^D]) in the limit 
of infinite energy. In this point, we can see a unifi- 
cation of quantum-physics description of state of ideal 
gas/radiation with the general-relativistic description of 
space-time filled in with a gas and radiation. 

Although we demonstrated that the efficiency of the 
URMS cooling is low when the condition \ 3> W is valid, 
a certain energy still escapes away and this energy loss 
obviously leads to shrinking of the URMS until its proper 
physical radius becomes Rb = Rg- For Rb = R g , no reg- 
ular energy loss is possible and the equilibrium set by the 



identity of Eqs. (|49|) and (|70|) is exact. This equilibrium 
seems to be the law of Nature, supposed by Eddington, 
which prevents very massive stars from the reductio ad 
absurdum. According to our result, the final stage of an 
object with a mass exceeding the Oppenheimer-Volkoff 
limit is the stable object consisting of a matter and ra- 
diation and having radius exactly equal to R g . In other 
words, the outer border of the object is situated just at 
the horizon of events. 

The concept of the space-time surrounding the URMS 
does not differ from the classical concept of the space- 
time surrounding the non-rotating, Schwarzschild black 
hole. Also the existence of the horizon of events is pre- 
dicted within both concepts. In our concept, this horizon 
cannot be, however, regarded as an absolute, insuperable 
barrier for any matter or radiation to escape. Just at the 
horizon, there are permanently situated the gas and ra- 
diation having the internal energy approaching infinity. 
In this context, it is necessary to note that the derived 
equations are related to a perfectly quiet fluid with no 
local phenomena. But in the reality, we do not know any 
perfectly quiet object. Also on the URMS surface, some 
fluctuations of quantities of state can be expected. Due 
to these fluctuations, some matter cannot be excluded 
to occur above the horizon. If this happens, an emission 
of a radiation can be expected. It must, of course, be 
extremely red-shifted, but its eventual detection cannot 
be excluded. In this aspect, our concept of the URMS 
is different to the concept of classical black holes. (We 
started to use the term " URMS" , instead of the tradi- 
tional "black hole", especially because of this reason.) 

To explain the basic principles of our new concept of 
the extremely hot compact object more transparently, 
we constrained our considerations to the simplest case of 
spherically symmetric, non-rotating object. Despite this 
simplicity, we believe that the concept of URMS will be 
helpful to explain some observational phenomena related 
to the concerning objects. 
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